Non-deterministic noiseless amplification via non-symplectic phase space 

transformations 
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We analyse the action of an ideal noiseless linear amplifier operator, g° ta , using the Wigner 
function phase space representation. In this setting we are able to clarify the gain g for which a 
physical output is produced when this operator is acted upon inputs other than coherent states. 
We derive compact closed form expressions for the action of N local amplifiers, with potentially 
different gains, on arbitrary iV-mode Gaussian states and provide several examples of the utility of 
this formalism for determining important quantities including amplification and the strength and 
purity of the distilled entanglement, and for optimising the use of the amplification in quantum 
information protocols. 



I. INTRODUCTION 

Optical quantum communication has resulted in nu- 
merous protocols that achieve classically impossible tasks 
including teleportation [J Q , quantum key distribution 
[H, H| and super-dense coding 0, @. Furthermore, sev- 
eral of these have seen experiments ranging in sophisti- 
cation from proof-of-principle demonstrations |7 H1 1 to 
implementations approaching real world conditions 1 2| — 
[l5| . One of the great challenges that stands between 
these schemes and the realisation of large scale quantum 
information networks is the necessity of preserving often 
fragile quantum states in the presence of losses and other 
decoherence. A device that allowed for amplification to 
combat such effects would be extremely useful, however 
the laws of quantum mechanics themselves conspire to 
enforce a noise penalty whenever such an operation is 
attempted 

An ingenious recent approach is to circumvent these 
limits by designing devices that achieve genuinely noise- 
less amplifications in a non-deterministic but heralded 
manner [l7j . This noiseless linear amplifier (NLA) has 
been the subject of considerable theoretical jLSl - l3l| and 
experimental j32Tr39| work. Applications in quantum key 
distribution (QKD) with both continuous variable (CV) 
[lH [I?], [28| and discrete variables fDV) l20ll23j have been 
considered as well as error correction [21[. In the litera- 
ture one finds two kinds of analysis. In the first place one 
can consider the ideal amplification operation, which is 
to implement g a ^ a . In the amplification regime (g > 1) 
this is an unbounded operator, however for any particu- 
lar input state one can always write down a new opera- 
tor of the form Il7v<? ata where LT^r is a projector onto the 
subspacc spanned by the first N + 1 energy eigenstates. 
This operation will result in an arbitrarily good approx- 
imation of an ideal NLA as N increases at the price of a 
decreased, but finite, success probability. In the second 
kind of analysis works thus far have utilised particular 



linear optics implementations such as that of the original 
proposal [1 711 or those based upon photon addition and 
subtraction]!! [U EJ . 

Here we will adopt the first approach, and focus on 
gaining a greater insight into the properties and applica- 
tions of the operation g a1a . In Section [XT] we will derive 
the action of the operation on an arbitrary state via the 
Moyal product and show that it has the unusual property 
of being a Gaussian but non-symplectic map. Further- 
more we address the scenarios in which, dependent upon 
the state to be amplified, the NLA fails to transform 
into a physical output in the limiting procedure described 
above. Some necessary concepts in Gaussian quantum in- 
formation are introduced in Section lull In Section HVl we 
apply our formalism to obtain compact analytic expres- 
sions for up to N amplifiers acting upon N-mode Gaus- 
sian states. We give specific examples for one and two 
mode cases and comment more rigorously on the physi- 
cality of the operation dependent upon the input state. 
In Section [V] we consider in more detail the distribution 
of EPR entanglement through a general Gaussian chan- 
nel which is the situation most relevant to continuous 
variable QKD. We show that previous attempts to repre- 
sent the action of the NLA and an effective channel of the 
same form but with different parameters are in general 
insufficient. In a surprising example we show that in the 
for large gains the NLA has the effect of transforming an 
attack on one half of an EPR pair into an attack upon 
the other half. Finally in Section WT\ we conclude. 



II. NLA AS A NON-SYMPLECTIC OPERATION 



We can define ideal linear amplification in terms of 
coherent states as 



\ga) (ga\ = T(\a) (a\). 
We can realise T by 
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where 



Tjv(p) =p 2 N 9 a a n NP Tl N g a 



(3) 



where the constant p^ is chosen to make the operation 
physical and II jy is the projection operator defined ear- 
lier. For most of this work we will ignore p^, it is impor- 
tant to say a few words about it at this point. To achieve 
ideal linear amplification over the entire harmonic oscil- 
lator Hilbert space requires lim.7v-i.00 Pn = as the g a ^ a 
operator has an unbounded spectrum. However, in any 
realistic experiment there will be some bounds within 
which it is assumed that the experiment is being per- 
formed. First there will be assumed some energy bound 
which can be thought of as a truncation of the Hilbert 
space. This results in accepting a non-unit fidelity with 
the theoretically ideal amplifier for states which have a 
component outside this bound. The choice of truncation 
is somewhat arbitrary, but will generally be determined 
by the energy limits of the experiment. We can think of 
a sequence of operations which is indexed by the largest 
energy eigenstate which is allowed, which we will call N. 
For any finite N there is a finite non-zero pm which one 
can choose for Yjy. As N grows, ppf must reduce. In the 
limit as N approaches infinity, we recover the ideal oper- 
ation and Pn tends to zero. The prediction we make by 
ignoring p^ in the theory is the state resulting from this 
limiting case. However we emphasise that the existing ex- 
periments have already shown that for low energy input 
states approaching the limiting case of ideal operation is 
achievable without prohibitively low success probability. 
The g a a operator in Wigner space is 
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where we have chosen h = 2 and used the identity a) a = 
j(x 2 + p 2 — 2). Whilst we are interested in the cases 
where g > 1, there is no such restriction needed for the 
calculations we will perform here. 

To act the Wigner representation of the amplifier on 
an arbitrary input state requires that the operator prod- 
uct be performed in the Wigner representation. This is 
achieved by way of the Moyal Product [4{| which we will 
denote *. If we take the input state p whose Wigner 
function is W p then the Wigner function for the output 

amplified state p' = g°^ a pg°^ a is, 



W. 
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We will later show that writing the action of the NLA 
in this form allows the calculation of compact analytic 
results for the class of Gaussian states but already the 
Wigner representation sheds some light upon the unusual 
properties of this operator. 

An interesting point noted in the original proposal fl7j 
is that the NLA will not produce a physical output when 



acted upon certain input states with certain gains. This 
does not come as a complete surprise given the state de- 
pendent manner in which the transformation is defined, 
namely over the coherent states. The Wigner representa- 
tion allows us to understand both the question of physical 
convergence and Gaussianity of the transformation. 

The expression Eq|4]is clearly a Gaussian operation in 
that it is an exponential quadratic in the phase space 
variables however for g > 1 the expression is a convex, 
unbounded function in phase space. Nonetheless when 
the NLA is acted upon a particular input state and the 
expressions are combined under the appropriate phase 
space convolution, EqH] tells us that if the input state 
is sufficiently 'small' (i.e. in terms of the decay of it's 
Wigner function in phase space) relative to the gain of 
the NLA then the overall output will have a concave, nor- 
malisable phase space distribution. That is, the limiting 
state linijv-Mnf T n (p) is a well defined physical state for 
this particular situation. Furthermore in this case the 
NLA will preserve the Gaussianity of the input state in 
the infinite limit. In section IIVI we will make rigourous 
this argument for the restricted class of Gaussian states, 
give exact criteria for the convergence of the output for 
a given input. 

Finally it is straightforward to see that the opera- 
tion although mapping Gaussian input to Gaussian out- 
puts, it is not symplectic. This is not in violation of 
the well known Stone- von Neumann theorem however as 
the operation also fails to be unitary g a a = (g a => 
giJa^gcJay 1 ^ c now t urn t the problem of evaluat- 
ing the action of the NLA in this phase space represen- 
tation upon the Gaussian states. First we will introduce 
some crucial properties and quantities of these states. 



III. SOME PRELIMINARIES OF GAUSSIAN 
QUANTUM INFORMATION 

Within continuous variable quantum information a 
great deal of attention is devoted to states with a Gaus- 
sian Wig ner function and operations that preserve this 
form 4l( . These so-called Gaussian states and operations 



are significant experimentally as they can be efficiently 
created and implemented and also lend themselves to an 
elegant theoretical description. This theoretical conve- 
nience comes from the fact that although such states live 
in in infinite dimensional Hilbert space they can be com- 
pletely characterised by their first and second moments. 

In particular if we start A^-modc Gaussian state, p, liv- 
ing in a tensor product of N infinite dimensional Hilbert 
spaces equipped with bosonic creation and annihilation 
operators Si, a\, ...djv, ajy we can use the corresponding 
quadrature operators £; = fij + a\ and Pi = i(a\ — a,) as 
phase space coordinates and completely characterise the 
state as follows: grouping the quadratures together in a 
vector r := (xi,pi, xn,Pn) we define any Gaussian 
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state by a displacement (mean) vector, 

d = tr(pr) 
and covariance matrix 

Ey =tr(p{(r i -d i ),(r 3 --d J -)}- 



where 



(5) 



(6) 



where {}+ is the anti-commutator. Matrices will be de- 
noted by boldface throughout and products of such terms 
should be interpreted as matrix multiplication. For a 
given square matrix to be a legitimate covariance matrix 
(CM) it must satisfy the uncertainty principle which in 
this formalism is the following positive semi-definitcness 
condition 143. 



where 
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which is called the symplectic form. 

The purity of a Gaussian state is obtained simply via 
the determinant of the CM 



H := tv(p 2 G ) 



1 



(9) 



Thus a Gaussian state is pure if and only if we have 
det(E) = 1. 

Likewise the entanglement of a Gaussian state is com- 
pletely determined by the second moments, in particular 
their symplectic spectra which we obtain by making use 
of Williamson's theorem [43j|. For any N-mode CM E 
there exists a symplectic diagonalisation given by 



N 



E = SWS T , W = 0A fe l 



(10) 



k=l 



where S is a symplectic matrix and Afc > 1 are the N 
symplectic eigenvalues of E. They can be computed ei- 
ther by solving an N th order polynomial where the co- 
efficients are so-called symplectic invariants [Hj or by 
finding the standard eigenspectrum of |iOS| where the 
absolute value is used in the operatorial sense Us- 
ing the fact that for symplectic matrices det(S) = 1 the 
physicality condition for a CM can be rewritten as, 



E > 0, A_ > 1 



where A_ is the smallest symplectic eigenvalue 

A common bipartite entanglement criteria is to con- 
sider the partial transpose (PT ) o f the density matrix 
with respect to one subsystem |45l . Hfjj, with entangle- 
ment (separability) corresponding to the non-physicality 
(physicality) of the resultant operator. For Gaussian 
states a similar results hold based upon the CM con- 
dition, namely that if we consider the partial transpose 
of a bipartite N x M mode CM 



E = (1a 95 T B )E(I, 



(12) 



(13) 



then entanglement corresponds to the non-physicality 
of E. It is straightforward to show that we always 
have E > Othcrcfore by Eq[TT] the criterion boils down 
to checking 



A_ > 1. 



(14) 



It should be noted that while Eq[TJ] is generally a nec- 
essary condition for bipartite separability it is only nec- 
essary and sufficient for 1 x N 14 71 and the bisymmetric 
class of N x M Gaussian states [48( . For such states how- 
ever one can quan tify the entanglement by means of the 
log negati vity [49( which can be computed for continuous 
variables [50( from the smallest PT symplectic eigenvalue 



£ = max{0, log A_} 



(15) 



Finally a large class of relevant operations and channels 
(squeezing, passive mode mixing, amplification and white 
noise channels) fall in the class of Gaussian operations. 
These can then be compactly implemented at the level 
of covariance matrix transformations. In this formalism, 
Gaussian unitary operations correspond to transforma- 
tions of the form E out = SSj n S T where S is a symplec- 
tic matrix. The NLA does not take this symplectic form, 
nonetheless we will derive analytic input output relations 
for the amplification of arbitrary Gaussian states. 



IV. AMPLIFICATION OF GAUSSIAN STATES 

It is inconvenient to evaluate the ★ operation in the 
x and p variables. It is more useful here to write this 
product down in terms of the characteristic functions 
(x) which are the Fourier transform of the corresponding 
Wigncr functions. Here we choose those variables to be a 
and b corresponding to transformed x and p respectively. 
With our choice of H = 2 (such that the vacuum noise is 
unity) the product in the Fourier transformed space for 
operators A and B is: 



(11) X AB(a,b) 



1 



(2tt)' 



da'db' X A{a-a', b-b')e l&Trta ' XB {a\ U) 



(16) 

where a = (a, b) T , a! = (a',b') T . This formalism carries 
over to the multidimensional case fairly naturally with 
the matrix 1~2 extended using the direct sum. The char- 
acteristic function for multidimensional Gaussians is 



XG 



(17) 



The Moyal product of two Gaussian will then result in 
another Gaussian. Thus we can simply read off the out- 
put covariance matrix and hence overall state from the 



4 



output characteristic function. Also, we can consider act- 
ing amplifiers on all modes and the amplification trans- 
formation extends naturally from Eq|4] where we rewrite 
as 



G w (x,p) =cxp{G-V+p 2 )} 
we have defined the matrix 

G = 




(18) 



(19) 



and extend G to the multimode case by taking the direct 
sum of all the individual modes. We can take the limit 
as g — > 1 for modes which have no amplification. This 
limit constitutes a delta function as expected. 

To conjugate the density operator by the action of an 
ideal amplifier the Moyal product must be applied twice. 
For a Gaussian density operator with mean vector d and 
covariance matrix E this computation gives two equa- 
tions for the mean vector and covariance matrix, one from 
each application of the product. 

To compute the relationships for the mean vector and 
covariance matrix we will utilize the result of the gaussian 
integral 



dx\dx2 ■ ■ ■ dx n e 



"A-x+b-x K 



(20) 



where this equation holds for matrix A and vector b 
possibly complex (note that the transpose is taken and 
not the adjoint). We can ignore the proportionality factor 
here as it can be considered part of the pn terms we 
defined in ScctionTIl Since we are only interested in the 
convergent state this factor is not relevant. 

We will start by computing the left product by the 
g a a operator p\ = g a a p (which is not a physical state) 
using the Wigncr representation and the Moyal product. 
Substituting in the definitions for the operators results in 
the integral (ignoring constant proportionality factors) 



yy^ _ J ^ a ' e (a-a') r G(a-a')+ja' T f2a+id T a'-ia' T £a' 

(21) 

We can now rearrange the polynomial to obtain an ex- 
pressions more like that of the standard Gaussian inte- 
gral. 

Wpi = e 5 aTGa J da ' e «d T +a T n)-a T G)a'-ia' T (E-G)a' 

(22) 

Upon evaluating this integral we find that the covariance 
matrix transforms as 

Si = (in — G)(E — G) _1 (in + G) - G (23) 

and the mean vector transforms as s 

di = (in - G)(E - G) _1 d (24) 



Following a similar calculation for right multiplying 
by the amplification operator transforms the covariance 
matrix as 

S OU i = (in + Ei)(£i - G)-\in + Si) + Si (25) 



and the mean vector transforms as 



di - (i« + Si) (Ei -G) _1 di. 



(26) 



In the single mode case, G will be proportional to the 
identity and hence this can be used to greatly simplify 
the expressions, but we will work with the case of G not 
being proportional to the identity and hence will include 
multimode cases where different gain variables may be 
used on different modes. 

To simplify these expressions it is important to note 
that 



(in-G)(in + G) =I-G 2 



(27) 



This holds as G is diagonal and for each local block the 
two entries are the same. This gives f]G = Gfl and 
hence the above relation. 

We will now dissect the terms in S ollt and evaluate 
them: 

(Ei-G)- 1 = (itt + G)- 1 

[(£ - Gr 1 + (in + g)- 1 - (in - g)- 1 ] _1 

(in-G)" 1 (28) 

in + Si = (in - g) [(s - g)- 1 + (in + g)- 1 ] 

(in + G) (29) 

in -Si = (in - g) [-(e - G)- 1 + (in - g)- 1 ] 

(in + G) (30) 
Substituting these results back gives 
E otlt = [(S - G)- 1 - 2(G" 1 - G)- 1 ] 1 G (31) 
This can be further rearranged to give 

E OMt = G- 1 (G- 1 + G-2S)" 1 (S-G) 



G(G- 1 + G-2S) (E-G" 



(32) 



If instead of using the matrix G with diagonal elements 
of the form we use a matrix 



9 
9 



(33) 



extended to many modes also using the direct sum, then 
this matrix equation becomes 

S out = g [g 2 + 1 - S(g 2 - l)]" 1 

[S(g 2 + 1) - (g 2 - 1)] g- 1 (34) 
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This equation can be further simplified if we replace the 
matrix of linear gains g to a matrix of logarithmic gains 

> = ('o 9 i„° 9 ) <»> 

Y, out = (coshl - £sinhl) _1 (£ coshl - sinhl). (36) 

Using the same relationships found above, we can now 
substitute into the expression for d out and get 

d out = [2(E - G)(G - G" 1 )" 1 -l] _1 d (37) 
and using the g matrix form gives 

d out = 2g [g 2 + 1 - £(g 2 - 1)] _1 d (38) 

and using the logarithmic gain form gives 

d out = (coshl - £smhl) _1 d (39) 

We will now use the results from this method to illu- 
minate how we can easily compute the action of the NLA 
on Gaussian states in one and two-modes. 

A. Single-mode states 

In general one can unitarily transform one mode Gaus- 
sian to remove cross-correlations between the quadra- 
tures and hence write CM's of diagonal form, 

with relevant examples including thermal, squeezed and 
coherent states. Considering such states with input mean 
vector d = ((x) , (p)) T we transform via [38] to find 

/ 2g(x) \ 

«W = (41) 

V V p + l-gi(V p -l) ) 

Naturally for coherent states (V x = V p = 1) the above 
expression reduces to d out = gd but for thermal states 
(V x = V p = V > 1) the amplification is no longer linear 
in g. Instead it increases rapidly and becomes infinite at 
a maximum value of the gain before becoming negative 
and as we will see unphysical. 

While both thermal and coherent states are phase 
symmetric, this is not true of squeezed states (V x = 
V~ l , V p = V > 1) and thus the two quadrature displace- 
ments amplify differently. While the anti-squeezed dis- 
placement transforms identically to the thermal case the 
squeezed quadrature displacement amplifies sub-linearly 
in g. These results are plotted in FigQ] 

Turning to the output covariance matrix we can sub- 
stitute int Eq|34]to obtain, 

/ T4 + l+g 2 (K-l) n \ 

^ NLA ~ n v P +i+g 2 (v P -D ■ ^ 2 > 

\ U V p +l~g 2 (V p -l) J 



2.4 
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FIG. 1: (colour online) Expectation values of output quadra- 
tures X — {x,p} for input coherent (black, dashed), thermal 
(blue, solid) and squeezed (red, x dot-dash, p crosses) states. 
For all plots V = 1.5 and d = (1, 1). 



Following a similar pattern the thermal state variance 
and that of the anti-squeezed quadrature for the squeezed 
state become infinite at the same gain for which the 
squeezed variance vanishes. In other words a thermal 
state thermalises further whereas a squeezed state is 
squeezed further FigEl Note that this is true regardless 
of the angle of squeezing. This pha se insensitive squeez- 
ing property was addressed in [26| where the potential 
causal paradoxes were resolved. 




1 1.2 1.4 1.6 1.8 2 



g 

FIG. 2: (colour online) Quadrature variances Vx, X = {x,p} 
for input thermal (blue, solid) and squeezed (red, x dot-dash, 
p crosses) states. For all plots V = 1.5 and d = (1, 1). 

This behaviour is in accord with our intuition from the 
first section in which we saw evidence for a critical gain 
at which the output state was described by a flat phase 
space distribution (infinite variance) and beyond which 
was unphysical (negative variance). 
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By considering the output on the single mode case we 
can derive the necessary and sufficient condition for con- 
vergence to a physical output state. Solving for the sin- 
gularity in the output CM Eq02]we find that for a single 
NLA applied to one mode of a multi-mode Gaussian state 
the maximum physical gain is 



V + l 
V-l' 



(43) 



where V is the variance of the input to the amplifier. 
Note that in the multi-mode case where more than one 
NLA is present this condition will be necessary but no 
longer sufficient to guarantee a convergent output. 



B. Two- mode states 

The canonical bipartite Gaussian state is the two-mode 
squeezed vacuum or EPR state which, among other ap- 
plications, is the building block for quantum tclcporta- 
tion and the theoretical analysis of continuous variable 
quantum key distribution. In the number basis it has 
the form, 



\EPR) = V^X 2 



(44) 



where x ranges from for an unsqueezed vacuum up 
to unity for an infinitely squeezed, maximally entangled 
state. 

The effect of the NLA in this context has already been 
considered in [21] where distillation in the presence of 
loss was demonstrated and it has also been shown to 
benefit key distribution over general Gaussian channels 
pM [271 . 12H ] . Here we will analytically re-derive the pre- 
vious results but with much less effort via our new for- 
malism. 

An arbitrary Gaussian channel can be described by a 
transmission T and a thermal noise parameter, for in- 
stance the variance of the enivronment Ve- 

The covariance matrix of an EPR state with mode B 
distributed through such a channel is given by, 



Sin — 



V A h 

CAB O z 



CAB O z 
V B h 



(45) 



with a. 



[1,0; 0,-1] and 

1 + x 2 



V A 
V B 
cab 



l-X 2 



= T- 



.1 



X 



l-X 2 



(1 - T)V E 



X 



1 



(46) 



Again applying Eg IM1 where an NLA of gain g is applied 
to mode B we find an output covariance matrix of the 



same form as Eq|45]but with entries related to the inputs 
via 

V'a = ^ [V A (V B + 1)-T(V 2 -1) 

-1 + VI)-V a {-1 + Vb))} 



where 





+ 9 2 




1 




N 




1 


c'ab 


N 




N 



[V B + l + g 2 {V B -1)] 



2 9x /T(V 2 - 1) 



l-g 2 (V B -l) 



(47) 



(48) 



which must be positive and non-zero and hence gives the 
constraint on g outlined in the previous section. 

We can now analyse the output states in terms of their 
entanglement and purity. We know from(l7l] that for an 
initially pure EPR state, as the NLA gain reaches it's 
maximum value the resultant output will tend towards a 
pure maximally entangled EPR state. The same results 
were also shown for pure loss channels where, crucially, 
no entanglement is generated between the transmitted 
mode and the environment. As soon as there is decoher- 
ence any amount of decoherence the NLA wtill start to 
distill correlations between the amplified mode and the 
noisy environment as well as the other arm of the EPR. 
Thus it becomes impossible to distill maximal entangle- 
ment within the allowable gain range. We will consider 
the entanglement distillation in the presence of Gaus- 
sian decoherence via the logarithmic negativity given in 
EqlTSl For two-mode states with CM of the form Eql45l 
the symplectic eigenvalues of the PT state are given by, 



\ 



A± J A 2 -4det(S) 



(49) 



where A := a 2 + b 2 + 2c 2 and det(S) = det(S) = (ab - 
c 2 ) 2 . In Fig [3] we plot the logarithmic negativity as a 
function of gain for the ideal channel and for increasing 
levels of decoherence. 

We paramaterise the channel by fixing the initial vari- 
ance (Va) (or alternatively the EPR paramater x) 01 
the EPR to be distributed and the variance of a ther- 
mal environment mode Ve and then interacting the two 
modes on a beamsplitter of varying transmission. A per- 
fect channel corresponds to T = 1 with an increasingly 
mixed output for smaller values. For the cases plotted we 
have Ve < Va so as the transmission decreases so does 
Bob's variance leading to a higher maximum allowable 
gain. These results show the entanglement increasing 
monotonically as a function of g. While a perfect chan- 
nel allows for maximum distillation, modest amounts of 
loss and noise swiftly decrease the maximum distillable 
entanglement and the rate of increase with g. 

We are also interested in the purity of our final state, 
and by considering Eq(5] it is swiftly apparent that the 
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g 



FIG. 3: (colour online) Logarithmic negativity of an amplified 
EPR state as a function of gain for varying levels of channel 
decoherence. We paramaterise this by fixing an initial EPR 
strength of x = 0-4 (Va ~ 1.4) and a thermal environment 
of variance Ve = 1.1 and mixing the two on beamsplitter 
of decreasing transmission T. Plotted curves are for T = 1 
(blue), T = 0.8 (red) and T = 0.5. 

purity decreases monotonically with gain and that 

lim n = 0. (50) 

g^gmax 

We plot the purity of the output states corresponding to 
those shown in Fig[3] demonstrating the degradation of 
purity as a function of gain. We see that the decay is 
initially gentle and then rapidly decreases as the maxi- 
mum gain is approached. Thus we are left with a further 
restriction upon the maximum distillable entanglement 
if we would also like to maintain high levels of purity. 




g 



1.5 2.0 2.5 3.0 



FIG. 4: (colour online) Purity of an amplified EPR state as a 
function of gain for the same input parameters as Fig(3] Once 
again Plotted curves are for T = 1 (blue), T = 0.8 (red) and 
T = 0.5 

These results demonstrate the competing considera- 
tions of entanglement strength and purity that must be 
taken into account when choosing the NLA gain to be 
applied. In the previous discussion ignored the final de- 
gree of freedom available to Alice and Bob, namely the 



strength of the input EPR. For example in a situation 
where the purity of entanglement is of great importance 
Alice and Bob can attempt to improve their protocol by 
starting with weaker entanglement and amplifying it fur- 
ther after the channel. As a further example of the utility 
of our formalism is the ease with which one can opti- 
mise over input parameters to maximise a desired op- 
erational quantity. As a demonstration we will plot the 
maximum achievable fidelity between the de-cohered and 
subsequently amplified EPR state and a pure target EPR 
state of a certain strength. We plot the results as a func- 
tion of the target EPR parameter \t where both input 
EPR strength and the NLA gain have been optimised 
over. The fidelity between two-mode Gaussian states of 
zero mean is given by [5l| , 




FIG. 5: (colour online) Maximum fidelity between a dis- 
tributed EPR state after amplification and a pure target EPR 
state as a function of the target EPR strength \t for varying 
levels of channel decoherence. We paramaterise the channel 
by fixing a thermal environment of variance Ve = 1.01 and 
mixing the two on beamsplitter of decreasing transmission 
T. Plotted curves are for T = 1 (blue), T = 0.9 (green) 
and T = 0.5 red). For all points the NLA gain and input 
EPR strength are simultaneously optimised over. To illus- 
trate the improvement due to amplfication we also plot the 
performance in the absence of an NLA (dashed lines) for com- 
parison. 
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where the input covariance matrices are in standard 
form with entries Va, Vg, cab and V' Al Vg, c' AB respec- 
tively. This expression differs slightly from that given 
in [5l| as we have a different noise convention with our 
vacuum normalised to 1 instead \. 



V. EFFECTIVE CIRCUIT 

In the previous section we derived the effect of the 
NLA upon an EPR state transmitted through an arbi- 
trary Gaussian channel and calculated the strength and 
purity of the resultant ent ang lement. The same situation 
was also considered in [25|,[27j[28| f° r the purposes of Con- 
tinuous Variable Quantum Key Distribution (CVQKD) 
where the authors paramaterised the results in terms of 
an effective combination of a different EPR state and 
channel. However if we pursue this representation in de- 
tail we find that although such a description is helpful as 
a tool for conceptualising how the NLA will affect com- 
munication rates it is insufficient to describe the true 
nature of correlations with an eavesdropper (Eve) . 

One can attempt to solve for a set of effective param- 
eters T',£',x\ where the £ is another noise parameter 
more commonly used the CVQKD literature. It is de- 
fined by setting the total variance added by the channel 
environment to be Ve = 1 ~- l T 1p T ^ . From this formula it 



1-T 

is apparent that the added channel noise has been split 
up into a component due to loss and the so-called excess 
noise £. In these papers the output of the NLA corre- 
sponds to a scenario where Alice created a stronger EPR 
initially and transmitted it through an effective channel. 
Setting the equations in Eqs l46l and 1471 equal we obtain, 



A 

r 



T + 5 (1 ~_ g2)) v 



4.g 2 T 



(-2 + (-1 + ,g 2 )T(-2 + 0)(-2 + (-1 + g 2 )TO 



--(-2 + (-I + .g 2 )T(-2 + £))£. 



(53) 



The secret key rates of the previous works are of course 
correct as they only depend upon the reduced covariance 
matrix shared by Alice and Bob, however this interpre- 
tation is not always valid and it turns out does not fully 
capture the correlations generated by the NLA. 

In general one assumes that all of the observed noise 
originates from the eavesdropper interactions, in other 
words that Eve purifies the state. For a noisy Gaussian 
channel where Bob's mode is effectively mixed with a 
thermal state a valid purification is for Eve to have her 
own EPR state and interact one arm with Bob's mode. 

Consider the symmetrical case where Alice and Eve 
create identical EPR pairs and interact them upon a 
50:50 beamsplitter. If Bob applies an NLA and we ask 
about his correlations with Alice and then Eve then the 
previous analysis results in a contradiction. If we con- 
sider the Alice-Bob channel we we see the effective trans- 



mission increase, but if we consider the Eve-Bob channel 
exactly the same result should hold. Clearly the beam- 
splitter ratio cannot simultaneously increase and decrease 
and we arrive at a contradiction. 




FIG. 6: Effective channel transmission as a function of gain, 
as given by Eq |53l Input parameters are the same as Fig[3] 
with T = 0.5. The curve is only plotted for gains less than 
the maximum allowed value but the effective transmission still 
exceeds the maximum sensible value for a beamsplitter. 

In fact if we plot the effective transmission given in 
Eql53l as a function of gain, Fig|6l we see that for large 
but allowable gains they cease to make physical sense. In 
particular the effective channel transmission can surpass 
unity, indicating that the NLA can not always effectively 
be equated with a beamsplitter interaction. Note that al- 
though this effective paramaterisation has broken down, 
key rates calculated by considering only the reduced co- 
variance matrix are still valid as the entropies are inde- 
pendent of Eve's particular purification. Nonetheless we 
are interested in ascertaining the exact form of the inter- 
action given by the NLA. 

The ease with which our method can be adapted to a 
multi-mode picture allows us to straightforwardly answer 
this question by explicitly including Eve's modes in the 
calculation and explicitly analysing the correlations. We 
consider the situation where Eve makes an entangling 
doner attack as per the upper panel of FiglZl mixing one 
arm of her own EPR with Bob's mode on a beamsplitter 
of transmission T. Bob subsequently applies an NLA 
before detection. The initial 4-mode covariance matrix 
looks like, 



( V l 2 CAB O z 

CAB <7z V I 2 



V o o 






V E Ha 

CE1E2 &z 






ceiei cr z 
V E h 



(54) 



with cab = VV 2 — 1 and ceiei = \JV^ — I. The final 
CM is obtained by enacting a beamsplitter transforma- 
tion between Bob's mode (B) and the first eavesdropping 
mode (EI) that is, £ BS BB1 (T)£BS^ E1 (T) and then 
substituting this into Eq |34l This results in, 
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FIG. 7: Equivalent eavesdropper attack in the presence on 
an NLA: a) In the original scenario Alice and Eve each cre- 
ate EPR which they mix on a beamsplitter of transmission 
T with the transmitted mode being sent to Bob, who uses 
an NLA. b) The correlations generated by this are identical 
to an equivalent scenario where Eve mixes her unused EPR 
arm with Alice's mode on a beamsplitter of transmission Ta 
while the original beamsplitter changes transmission to a dif- 
ferent value Tb- Both effective EPR pairs also increase in 
entanglement. 



'NLA 



where, 



/ V A la 

CAB &z 
CAE! &z 
\ CAE2 I2 



cab 0"z cae\ cr z 

Vb h c B Ei h 

cbei h Vei I2 

CBEI Oz CE1E2 O z 




V A 
V b 
Vei 

V E 2 



= V + T+(1-T)VV E 
+ g 2 {V-T-{l-T)VV E ) 
= TV + (1 - T)V E + 1 
+ g 2 (TV + (1-T)Vb)-1) 
= (1 — T)V + TVe + We 
+ g 2 {{l-T)V + TV E -VV E ) 
= V E + l + T(VV E -l) 
+ g 2 (V E - 1 - T(VV E - 1)) 



c AB = 2g^T(V 2 - 1) 
caei = -(V E + lW(l-T((V 2 -l) 
+ g 2 (V E -l-T(VV E -l)) 

cbei = 2gy/(l-T)T(y-V E ) 

2g\ 



CBE2 



CE1E2 



CAE2 



-T)(y|-i) 

2g ] f(l-T)(V 2 -l) 



(g 2 1)^(1 -T)T(V 2 -1)(V 2 -1) (56) 



however if we consider the correlations with the eaves- 
dropper we notice an extremely unusual feature. Ordi- 
narily for any Gaussian channel acting solely upon Bob's 
side the correlation between Alice's mode and the eaves- 
droppers non- interacted mode, cae2, is identically zero. 
However examining this term in Eq l56l we find that when- 
ever g ^ 1 there are correlations despite the fact these 
modes never interacted. 

Such correlations could never be reproduced by an ef- 
fective setup of the form of the top panel of Figj7]so we 
are motivated to construct an equivalent setup including 
an eavesdropping attack on both Alice's and Bob's modes 
as per the bottom panel of Fig|7] This would produce a 
CM given by, 



BSae2(T a )BSbei(T b )'Z 

BS B£l( r -B) BS AE2(T.A) 



(57) 



If we equate EqJSH] and EqJST] some lengthy algebra does 
result in a unique solution for the effective scenario where 
Eve interacts one EPR mode with Alice and one with 
Bob. The expressions for these effective parameters are 
given in detail in Appendix A. 




The elements of the reduced CM shared between Alice 
and Bob arc exactly the same as in the previous section, 



FIG. 8: Effective channel transmission Alice's (Ta) and Bob's 
(Tb) side as a function of the NLA gain. The input param- 
eters are T = 0.5 and Ve = 1.1. As the NLA gain increases 
the channel on Alice's side worsens (7a decreases) whereas 
the channel on Bob's side improves (Tb increases) asymptot- 
ing to a perfect channel on Bob's side with the attack now 
exclusively on Alice's side. 

Considering the same initial state and channel as the 
previous section the effective channel parameters reveal 
an intriguing conclusion about the effect of the NLA. As 
the NLA gain increases the channel on Bob's side im- 
proves (i.e. the effective transmission increases) whereas 
the converse is true of Alice's side as shown in Fig|51 In 
fact for these parameters the attack on Bob's side disap- 
pears almost entirely. 

Finally we also see that both Alice and Eve's effective 
initial entanglement is increased though not equitably. 
The variance characterising Eve's effective EPR increases 
modestly as the gain approaches it's maximum whereas 
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FIG. 9: Effective entanglement for Alice (bottom panel) and 
Eve (top panel) paramaterised by their variance as a function 
of the NLA gain. The input parameters are T = 0.5 and Ve = 
1.1. The effective entanglement of both parties increases as a 
function of the gain. 



regimes resulting in physical outputs for input states 
other than coherent states, explaining this in terms of 
the relative divergence of the NLA and the target state 
in phase space. For Gaussian states we have derived 
compact analytic formula's for the action of up to N 
amplifiers upon A^-mode states. We have given explicit 
results for important examples in one and two modes 
and analysed the strength and purity of entanglement 
of EPR states through general Gaussian decoherence, 
demonstrating that our methods allow for the swift iden- 
tification of the best strategy for a given protocol. Finally 
we uncovered some intriguing effects on correlations be- 
tween two, two-mode entangled states with the counter- 
intuitive result that under amplification an interaction is 
moved from one side to another. 

Future work of great practical importance will be to 
investigate various proposals for implementations of the 
NLA, especially the promising post-selection approaches 
of [23, [28| , to determine corresponding success probabil- 
ities which are critical for protocols where the rate is of 
great importance such as QKD and metrology. Finally 
other results enabled by this work will include revisit- 
ing various canonical quantum information tasks such as 
cloning and state discrimination. 

Acknowledgements. — The authors would like to thank 
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Alice's diverges as expected from the previous thermal 
state results. 



VI. CONCLUSIONS 



In conclusion we have considered the phase space rep- 
resentation of an ideal NLA, gaining insight into the 



References 



[1] C. H. Bennett, G. Brassard, C. Crepeau, R. Jozsa, 

A. Peres, and W. K. Wootters, Phys. Rev. Lett. 70, 1895 [8] 
(1993). 

[2] S. Braunstein and H. Kimble, Phys. Rev. Lett. 80, 869 [9] 
(1998). 

[3] C. H. Bennett and G. Brassard (Proceedings of Inter- [10] 
national Conference on Computers, Systems and Signal 
Processing, Bangalore, India, 1984). [11] 

[4] F. Grosshans and P. Grangier, Phys. Rev. Lett. 88, 
057902 (2002). 

[5] C. Bennett and S. Wiesner, Phys. Rev. Lett. 69, 2881 [12] 
(1992). 

[6] S. Braunstein and H. Kimble, Phys. Rev. A 61, 042302 

(2000). [13] 
[7] D. Bouwmeester, J.-W. Pan, K. Mattle, M. Eibl, H. We- 



infurter, and A. Zeilinger, Nature 390, 575 (1997). 

A. Furusawa, J. S0rensen, S. Braunstein, C. Fuchs, 

H. Kimble, and E. Polzik, Science 282, 706 (1998). 

C. Bennett, F. Bessette, G. Brassard, L. Salvail, and 

J. Smolin, Journal of Cryptology 5, 3 (1992). 

F. Grosshans, G. V. Assche, J. Wenger, R. Brouri, N. J. 

Cerf, and P. Grangier, Nature 421, 238 (2003). 

A. M. Lance, T. Symul, V. Sharma, C. Weedbrook, T. C. 

Ralph, and P. K. Lam, Phys. Rev. Lett. 95, 180503 

(2005). 

M. Sasaki, M. Fujiwara, H. Ishizuka, W. Klaus, 
K. Wakui, M. Takeoka, S. Miki, T. Yamashita, Z. Wang, 
A. Tanaka, et al., Opt. Express, OE 19, 10387 (5). 
M. Peev, C. Pacher, R. Alleaume, C. Barreiro, J. Bouda, 
W. Boxleitner, T. Debuisschert, E. Diamanti, M. Dianati, 



11 



[42] 

[43] 

[44] 
[45] 
[46] 

[47] 
[48] 

[49] 

[50] 
[51] 



J. F. Dynes, et al., New Journal of Physics 11, 075001 [41] 

(2009) . 

X.-S. Ma, T. Herbst, T. Scheidl, D. Wang, 
S. Kropatschek, W. Naylor, B. Wittmann, A. Mech, 
J. Kofler, E. Anisimova, et al., Nature (2012). 
J. Yin, J.-G. Ren, H. Lu, Y. Cao, H.-L. Yong, Y.-P. Wu, 
C. Liu, S.-K. Liao, F. Zhou, Y. Jiang, et al., Nature 488, 
185 (2012). 

C. M. Caves, Phys. Rev. D 26, 1817 (1982). 
T. C. Ralph and A. P. Lund, Quantum Communication 
Measurement and Computing Proceedings of 9th Interna- 
tional Conference p. 155 (2009). 

P. Marek and R. Filip, Physical Review A 81, 022302 

(2010) . 

J. Fiurasek, Physical Review A 80, 053822 (2009). 
N. Gisin, S. Pironio, and N. Sangouard, Phys. Rev. Lett. 
105, 70501 (2010). 

T. C. Ralph, Physical Review A 84, 022339 (2011). 
S.-Y. Lee, S.-W. Ji, H.-J. Kim, and H. Nha, Physical 
Review A 84, 012302 (2011). 

D. Pitkanen, X. Ma, R. Wickert, P. V. Loock, and 
N. Lutkenhaus, Physical Review A 84, 022325 (2011). 
J. Brask, N. Brunner, D. Cavalcanti, and A. Leverrier, 
Phys. Rev. A 85, 042116 (2012). 

R. Blandino, A. Leverrier, M. Barbieri, J. Etesse, 
P. Grangier, and R. Tualle-Brouri, Phys. Rev. A 86, 
012327 (2012). 

C. Gagatsos, E. Karpov, and N. J. Cerf, Phys. Rev. A 
86, 012324 (2012). 

J. Fiurasek and N. J. Cerf, arXiv (2012), 1205.6933vl. 
N. Walk, T. Symul, P. K. Lam, and T. C. Ralph, arXiv 
(2012), 1206.0936v2. 

F. Ferreyrol, N. Spagnolo, R. Blandino, M. Barbieri, and 
R. Tualle-Brouri, arXiv (2012), 1205.6195vl. 
H.-J. Kim, S.-Y. Lee, S.-W. Ji, and H. Nha, Physical 
Review A 85, 013839 (2012). 

C. Navarrete-Benlloch, R. Garcfa-Patron, J. Shapiro, and 
N. Cerf, Phys. Rev. A 86, 012328 (2012). 

G. Y. Xiang, T. C. Ralph, A. P. Lund, N. Walk, and 
G. J. Pryde, Nature Photonics 4, 316 (2010). 
F. Ferreyrol, M. Barbieri, R. Blandino, S. Fossier, 
R. Tualle-Brouri, and P. Grangier, Phys. Rev. Lett. 104, 
123603 (2010). 

A. Zavatta, J. Fiurasek, and M. Bellini, Nature Photon 
5, 52 (2011). 

M. A. Usuga, C. R. Miiller, C. Wittmann, P. Marek, 
R. Filip, C. Marquardt, G. Leuchs, and U. L. Andersen, 
Nat Phys 6, 767 (2010). 

F. Ferreyrol, R. Blandino, M. Barbieri, R. Tualle-Brouri, 
and P. Grangier, Phys. Rev. A 83, 063801 (2011). 
S. Kocsis, G. Y. Xiang, T. C. Ralph, and G. J. Pryde, 
arXiv (2012), 1208.5881vl. 

M. Micuda, I. Straka, M. Mikova, M. Dusek, N. J. Cerf, 
J. Fiurasek, and M. Jezek, arXiv (2012), 1206.2852vl. 
C. I. Osorio, N. Bruno, N. Sangouard, H. Zbinden, 
N. Gisin, and R. T. Thew, arXiv (2012), 1203.3396v2. 
J. Moyal, Mathematical Proceedings of the Cambridge 
Philosophical Society 45, 99 (1949). 



C. Weedbrook, S. Pirandola, R. Garcfa-Patron, N. Cerf, 
T. Ralph, J. Shapiro, and S. Lloyd, Rev. Mod. Phys. 84, 
621 (2012). 

R. Simon, N. Mukunda, and B. Dutta, Phys. Rev. A 49, 
1567 (1994). 

J. Williamson, American journal of mathematics 58, 141 
(1936). 

A. Serafmi, Phys. Rev. Lett. 96, 110402 (2006). 

A. Peres, Phys. Rev. Lett. 77, 1413 (1996). 

M. Horodecki, P. Horodecki, and R. Horodecki, Physics 

Letters A 223, 1 (1996). 

R. Werner and M. Wolf, Phys. Rev. Lett. 86, 3658 (2001). 
G. Adesso and F. Illuminati, J. Phys. A: Math. Theor. 
40, 7821 (2007). 

G. Vidal and R. Werner, Physical Review A 65, 032314 
(2002). 

R. Simon, Phys. Rev. Lett. 84, 2726 (2000). 

P. Marian and T. Marian, Phys. Rev. A 86, 022340 

(2012). 



Appendix A: 4-mode equivalent circuit 

Here we derive in detail the parameters for the effec- 
tive circuit shown in the bottom panel of Fig|7] We start 
with two EPR pairs belonging to Eve and Alice, para- 
materised by variances V',V^ respectively. To calculate 
the necessary output CM's we need the NLA transform 
given by Eq l34l and the 4-mode version of the beamsplit- 
ter transform. This induces correlations between the two 
target modes and acts as the identity upon the others. 
For example the 4-mode beamsplitter matrix acting upon 
modes 2 and 3 with transmission T is, 



BS 2 , 3 (T) = 



( h 








f 









VT" 

T 




T 



\ 



/ 



(Al) 



Note that the minus sign that appears on one of the cor- 
relation terms is essentially a choice of convention, or cor- 
responds to a choice to swap which mode enters which 
port of the beamsplitter. 

In this purified version of the protocol the transmission 
through the channel simply corresponds to mixing modes 
B (2) and and El (3) via[7Q] 



S c h = BSs i _ei(T)£BS|J S1 (T) 



(A2) 



and the substituting this into EqJM] which directly gives 
the terms in Eq[56] We now calculate the equivalent 
circuit CM given by EqJST] which has the block diagonal 
form of Eq l45l with 



J 



V A 



V^+T A {V -V' E ) 



12 



V B = V E +T B (V'-V E ) 

Vei = T B (V E -V') + V 

V E 2 = T A (V E -V') + V 

cab = -y/(l ~ t a)(1 ~ T B ){V E 2 - 1) + VTaTb(V' 2 - 1) 

caei = -yj(l - T A )T B (V E 2 - 1) - s/TaJi - T B )(V' 2 - 1) 

cbei = s/{l-T B )T B (y E -Y) 

cbe2 = \Jt a {1 - T B )(V E 2 - 1) + v/(-l + T A )T B (l - V n ) 

CE1E2 = ^T A T B (V E 2 - 1) - y/(l - T A )(l - T B )(V' 2 - 1) 

cab2 = VO- - Ta)T a (V - V E ) (A3) 



We now wish to set this expression equal to Eqj56] si- 
multaneously solve for V' , V E , T A , Tb- In practice one an 
proceed by considering only a few terms and then check- 
ing that the solution satisfies all terms. By equating the 



variance terms one can swiftly solve for three of the pa- 
rameters in terms of the fourth and the input variables 
obtaining, 







T b 




\v E ) 





{-l-g 2 )TV E + (l+g 2 + (-l+g 2 )(-l+T)V E )v' +V {-1+T-Ve+TV' +g 2 (-1+T+Ve-TV')) 
(-l-V)(l+Vk)+2(l+T(V -Ve) + Ve)V +g 2 ((-1+V)(-1 + Ve)-2(-1+T(V -Ve) + V e )V) 

(-l-V)(l+VE) + (l+T(V-VE)+V E )V'+g 2 (l- + V'+VE(-l + (-l+T)V')+V(-l + VE-TV')) 
-l-TV + (-l+T)V E +g 2 (-l+T(V-V E )+V E ) 



(A4) 



Considering, say, the cae2 correlation term we hnd a 
quadratic with the two solutions corresponding to a per- 
mutation of the roles of Alice and Eve. Considering any 
of the other correlation terms will yield only one con- 
sistent solution. One further subtlety is apparent upon 
considering the Caei term. In particular the expression 
EcilSBl is always positive whereas the corresponding term 



m 

Eq|A3]chan 

ges sign depending upon the relative mag- 
nitude of V and Ve- The resolution to this is that for 
situations where the original parameters are such that 
Ve > V we must change the phase convention for the 
beamsplitter on Alice's side. This corresponds to mov- 
ing the minus sign in Eq lAll to the other correlation term. 
The two solutions for the remaining parameter V are, 



V > V, 



E 



g 2 (-l+V)(-l + VE)-(l + V)(l + VE)-VB^+4:C 
2(-l-TV +{-l+T)V E +g 2 (-1+T(V -Ve)+V e )) ' 

V = { (A5) 

g 2 (-l+V)(~l + VE)-(l + V)(l+VE) + VB2+iC y v 
2(-l-TV +(-l+T)V E +g 2 (-l+T(V -Ve)+V e )) ' V ^ VE 



wher 



B = l + V + V E + VV E +g 2 {-l + V + VE-VV E ) 

C = (-1 - TV + (-1 + T)V E + g 2 (-l + T(V - V E ) + V E )) ((l - g 2 ) TV E + V (l - T + V E + g 2 (-l + T + V E ))) 

I 



Direct substitution then confirms that for all input pa- 



rameters we have a unique set of parameters that yield 
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an identical CM to that created by the NLA. 



